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We present a number of arguments in favor of the suggestion that the Marshall- 
Peierls sign rule survives the frustration in the square-lattice Heisenberg antifer- 
romagnet with frustrating next-nearest-neighbor (diagonal) bonds (J% — J 2 model) 
for relatively large values of the parameter J 2 / 'Ji- Both the spin- wave analysis and 



O 

the exact-diagonalization data concerning the weight of Marshall states support the 



above suggestion. 
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Many years ago Marshall proved, using a Lemma due to Peierls, the well-known 
theorem determining the amplitude phases of the set of Ising states building the 
ground-state wave function of the spin | Heisenberg antiferromagnet on bipartite 
lattices [1]. Later on, Lieb and Mattis generalized the theorem to arbitrary site spins 
and bipartite lattices without translational symmetries [2,3]. 

Assigned to the J\ — J 2 model 



H = Hi + aH 2 = ^ SjSj + a ^ SjSj , J\ = 1 a = J2/J1 = J: 



2 , 



(nn and nnn mean that the summations run over the nearest-neighbor and next- 
nearest-neighbor, diagonal, bonds on the square lattice, respectively), the theorem, 
applicable to the case a < 0, says: 

(i) The lowest eigenstate of the Hamiltonian (1) in each subspace determined by 
the eigenvalue M of the spin operator Sf otal reads 

*m = E4 m V>, 4 M) >o. (2) 

m 

Here the Ising states \m) are defined by 

\m) = {-l) SA ~ MA \m 1 ) ® \m 2 ) ® • • • ® \m N ) , (3) 

where \m,i), i = 1, • • • , N, are the eigenstates of the site spin operator Sf (— Sj < 
m < s^, S A = EigA s i, M A{B) = J2ieA(B) m i, M = M A + M B . The lattice consists 
of two equivalent sublattices A and B. s^ = s, % — 1, • • • , N , are the site spins. The 
summations in Eq.(2) are restricted by the condition Y<iLi m i — M, N being the 
number of sites. 

A direct result of the property c^ > 0, valid for each m from the basis set (3), 
is the nondegeneracy of ^m, since it is impossible to build up other ortho normal 
states without using negative amplitudes c^f\ 

(ii) The total spin of ^m is S = \M\. Since each M subspace contains only 
eigenstates with S > |M|, a result of (ii) is the relation 

E{S) <E{S+1) , S>0 , (4) 

with E(S) being the lowest energy eigenvalue belonging to S. It is worth noticing 
that the property (ii) is characteristic for a larger set of Hamiltonians (1), because 
the relation (i)=>- (ii) can always be proved. 



In what follows we address the model (1) with frustrating diagonal bonds, i.e., 
a > 0. According to Lieb-Mattis definition [3], if a > 0, the square lattice is not 
bipartite because it is not possible to find a positive constant c > such that the 
conditions Jaa, Jbb < c, Jab > c are satisfied for all bonds of the square lattice. 
The goal of the present Letter is to give arguments in favor of the suggestion that the 
Marshall-Peierls sign rule (i) can survive the frustration for relatively large positive 
parameters a. 

For references, let us firstly recall the important steps leading to the sign rule (i) 
[3]. We suppose that the ground state of (1) in the subspace M (with an eigenvalue 
Em) is writen in the form (2). The Schrodinger equation for the amplitudes c^Jf' 
reads 

T, A nmC ( f ] = (e m -E M )c^ , (5) 

n 

where 

A nm = B nm -aC nm = -(n\H* y \m) - a(n\H% y \m) , (6) 

£ m (£m > Em) is the eigenvalue of H z in the Ising state \m). 

Further, notice that the variational function $m = H m Ic^ll 771 )) producing the 
Schrodinger equation 

Y. A nm\cl M) \ = (e m -E M )\cM\ , (8) 

n 

should also have an eigenvalue Em- Let us suppose that the matrix elements A nm 
satisfy 

Then Eqs.(5,8) lead to 

n n 

which is possible if, and only if, c^ > 0. Finally, the chance of some c^ to vanish 
is excluded, for if it were so, Eqs.(5,10) would imply that c^ = for all m. These 
are the main points leading to the sign rule (i). 

Now, notice that for a > the above reasoning is not applicable since the 
inequality (9) will be violated for some matrix elements A nm . Nevertheless, for 



small enough frustration parameters a > 0, one can argue that the Marshall-Peierls 
sign rule survives the frustration. Indeed, on rather general variational grounds the 
amplitudes cffl are expected to be continuous functions of the parameter a. The 
only point where one may expect some discontinuity of the amplitudes is the special 
point a = where the system (1) acquires a new quality, namely, frustrating bonds. 
So, let us assume that for some m the amplitudes satisfy c^f\a = + ) < 0. Then 
in the limit a j + , (9) will imply 

lim A nm = B nm > , (11) 

showing that the inequality (9) is satisfied at a = + . Thus, the same reasoning 
as in the case a < implies that the above assumption is wrong, i.e., all the 
amplitudes should be positive, c^\a = + ) > 0. Finally, since all the possible 
Ising configurations take part in the ground state at a = + with finite positive 
amplitudes, it is clear that the same property conserves in some region of positive 
a, as well. 

Spin-wave analysis: The self-consistent spin-wave approach used below [4,5] is 
predominantly addressed to low- dimensional spin systems. The conventional spin- 
wave technique is supplemented with a condition for zero sublattice magnetization, 
thus fulfilling, by hand, Mermin- Wagner's theorem [6] at finite temperatures, or the 
same requirement on finite lattices. The J\ — J2 model has also been analyzed by 
this method in a number of recent papers (for a review, see, e.g., Ref.7 and references 
therein). Here we omit the details and directly present the spin- wave ansatz due to 
the theory mentioned above 



i) sw ~ exp (]T V a+6+ k ) \Neel) . (12) 

k 

Here \Neel) is the classical Neel state. The weight factors w^ are defined by w^ = 
Vk/uk , t>k and uy_ being the well-known Bogoliubov coefficients 

2vl = (i- v iy 1/2 -i , 2i£ = (i-»£)- 1/2 + i, (13) 

r] k = -— — — — , 7 k = -(cosk x + cos k y ) , I\ = cos k x cos k y . 

1 + fi — au [1 — 1 kj ^ 

(14) 
The Bose operators or, 6r come from the Dyson- Maleev transformation and live on 
A and B sublattices, respectively. The prime over sums means that k vectors run in 



the small Brillouin zone. \x is the Lagrange multiplier used to imply the condition 
for zero sublattice magnetization. The renormalization factor U, renormalizing the 
frustation parameter a in the last equation, is a result of the Hartree-Fock decoupling 
corresponding to the theory (for details, see Refs.8-10). Within the theory, U is a 
ratio of two short-range correlators determined by the self-consistent equations. It 
has been shown in a recent study [10] that the spin- wave ansatz (12) gives an exellent 
fit to the exact-diagonalization data for the relevant quantities of the model (1) up 
toasi 0.45 , provided the quasiclassical limit (large site spin s^ = s) for the factor 
U (see, e.g., Refs.11,12) 
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Lr _ l-0.102/2 S + O[(2 g ) 

l + 0.158/2s + O[(2s)~ 2 ] 

is used. U — 1 corresponds to the linear spin-wave approximation [13]. 
Now, let us rewrite (12) in the form 

ijj sw ~ exp (- J2 w ( r )^R+r) \Neel) , (16) 

R,r 

Re.4 

where the pairing function w(r) is defined by 

2 
w ( r ) = tt Jl'^k cos kr . (17) 

iV k 

The vector r connects sites from different sublattices. From the structure of the 

latter state (16), it is clear that the sign rule (i) breaks if, and only if, the pairing 

function w(r) changes its sign for some vector r connecting two spins living on 

different sublattices. For the 4x4 lattice this is just the vector r = x + 2y (and 

the vectors related by the lattice symmetry ). A numerical calculation demonstrates 

that u>(x+2y) changes sign at a point practically coinciding with the related N = oo 

case (aM = 0.323 for U = 1). In the thermodynamic limit N = oo, we have checked 

numerically that it is just the function w(x+2y) which firstly becomes negative when 

a increases. Unfortunately, we have not succeeded in finding any analytical proof 

for this observation. Finally, it is interesting to notice that the curve representing 

the exact weight of Marshall states vs a for the 4x4 lattice is surprisingly well 

reproduced from the finite-size spin-wave analysis based on the ansatz (16) [10,14]. 

Exact-diagonalization data: Here we present exact-diagonalization results for the 

weight of non-Marshall states (which do not follow the sign rule (i)) in the exact 



ground-state wave function of the J\ — J 2 model, Eq.(l), on 4 x 4 and 6x4 lattices 
in the interval < a < 0.52, Fig.l (A detailed and varied analysis of that matter 
for the S = state is presented in Ref.14). Here we study the ground-state wave 
functions corresponding to total-spin quantum numbers S = 0, 1, 2, and 3. It is 
found a relatively wide region where the non-Marshall states lack at all (for each of 
the studied states). For example, the exact upper bound au in the case N = 4 x 4, 
S = is «m = 0.28 (this number seems to be slightly size-dependent as seen 
in Fig.l, 6x4 lattice). It is worth noticing that even in the strongly frustrated 
region a ~ 0.52 the weight of the Ising states violating the rule (i) is very small as 
compared to the weight of states which fulfill the rule. Further, a well pronounced 
tendency is that in the states with larger spin S the weight of Marshall states is 
larger, too. In addition, a more detailed study also shows [14] that the exclusion 
of the non-Marshall states from the exact ground states does not have any drastic 
effect on the spin-spin correlators and the other relevant quantities characterizing 
the model (1). It is interesting to compare the above picture to the one produced 
by the J\ — J2 antiferromagnetic chain, Fig. 2. As it is clearly seen, the behavior is 
quite different: First, the range where the sign rule is exactly fulfilled is considerably 
reduced. Second, for states with larger spins the weight of non-Marshall states is 
larger, a tendency which is just the opposite to the one seen in the two-dimensional 
case. Precisely at a — 0.5 the exact ground state is known to be the spin-Peierls 
dimer state, so that the Marshall-Peierls sign rule is exactly fulfilled. 

Concluding remarks: In conclusion, a number of arguments has been given in 
favor of the suggestion that in the J\ — J 2 model the Marshall-Peierls sign rule 
survives the frustration in a wide region up to the limit a«. Evidently, at this 
point some of the Ising amplitudes in the ground state should vanish and further 
become negative. This new quality of the ground state, however, does not sufficiently 
influence the quantities characterizing the long-range order in the system, which is 
rather natural in view of the picture drawn above. It is worth noticing that the 
number of non-Marshall states just after the Marshall point ckm, seen in the exact 
results, is very small, which is in agreement with the spin-wave picture showing that 
the violation of the sign rule starts from Ising states with a small number of reversed 
(in respect to the classical Neel state) sublattice spins. Further, the weight of Ising 



states violating the rule (i) remains extremly small up to the strong frustration limit 
a ~ 0.5. This fact has some importance for the construction of variational wave 
functions, for which the rule (i) is applicable even in the strongly frustrated region. 
It is also clearly seen that the ground states with larger spins S are, in a sense, 
less destroyed exactly in that case (two-dimensional J\ — J2 model) when the ground 
state is, most probably, macroscopically Neel ordered. In addition, it is interesting 
to notice that the sharp increase in the weight of non-Marshall Ising states with 
larger spins, Fig.l, is observed close to the same point a ~ 0.5 where, according 
to the renormalized spin-wave ansatz (15,16) (see Ref.10), the two-sublattice long- 
range Neel order is completely desrtoyed. Therefore, there are some indications in 
favor of the suggestion that the phenomenon discussed in the present Letter may, in 
principal, be effectively used in the finite-size analysis concerning the mechanism of 
breaking of the global rotational symmetry of the Heisenberg Hamiltonians [15,16]. 
Clearly, further work in this direction is needed. For such studies, the frustrated 
J\ — J2 Heisenberg antiferromagnet seems to be quite an adequate model, for the 
whole process of destroying of the Neel order can be suitably studied by changing 
the continuous frustration parameter J2/J1. 
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Captions of figures 

Fig.l: The weight of the non-Marshall Ising states in the lowest-energy eigen states 
with total spins S = 0,1,2, and 3 vs the frustration parameter J2/J1 for the s = 1/2 
square-lattice Heisenberg antiferromagnet with frustrating diagonal bonds for 4x4 
and 6x4 lattices. 

Fig. 2: The weight of the non-Marshall Ising states in the lowest-energy eigen 
states with total spins S — 0,1, 2, and 3 vs the frustration parameter J2/J1 for the 
s = 1/2 Heisenberg antiferromagnetic chain with frustrating next-nearest neighbor 
bonds for N = 16 and iV = 24. 



